Chain-induced effects in the Faraday instability on ferrofluids 
in a horizontal magnetic field 
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The linear stability analysis of the Faraday instability on a viscous ferrofluid in a horizontal mag- 
netic field is performed. Strong dipole-dipole interactions lead to the formation of chains elongated 
in the field direction. The formation of chains results in a qualitative new behaviour of the fer- 
rofluid. This new behaviour is characterized by a neutral stability curve similar to that observed 
earlier for Maxwell viscoelastic liquids and causes a significant weakening of the energy dissipation 
at high frequencies. In the case of a ferrofluid with chains in a horizontal magnetic field, the effective 
viscosity is anisotropic and depends on the field strength as well as on the wave frequency. 

PACS numbers: 83.60.Bc, 83.60.Np, 83.60.Wc, 83.80.Hj, 47.20.Gv, 75.50.Mm 
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I. 



INTRODUCTION 



Magnetic fluids (or ferrofluids) are colloidal dispersions 
of single domain nanoparticles in a carrier liquid. The 
fascination of ferrofluids stems from the combination of a 
normal liquid behaviour with the sensitivity to magnetic 
fields. This enables the use of magnetic fields to con- 
trol the flow of the fluid, giving rise to a great variety of 
new phenomena and to numerous technical applications 1 . 
One of the peculiarities of ferrofluids is the mutual influ- 
ence of the microstructure and the rheological properties 
of the fluid. Due to strong interparticle interactions var- 
ious aggregates can be formed in a ferrofluid^*^. On 
the one hand, the formation of the aggregates changes 
the effective viscosity of the fluid. On the other hand, 
the motion of the fluid influences the structure of the 
aggregates. 

In an applied magnetic field chains containing several 
particles are the favored form of an aggregate and are 
formed by the dipolar interactions. There is a similarity 
between the chains of dipoles and the macromolecular 
chains in polymer solutions. In both systems a network 
of chains is coupled with a viscous carrier liquid. This 
viscous coupling implies, that the relaxation character of 
the chain dynamics leads to a viscoelastic behaviour of 
the solution. The viscoelastic behaviour is reflected in 
a dependence of the stress tensor on the history of the 
system. The theoretical treatment of the viscoelasticity 
of polymer solutions is usually based on a phenomcno- 
logical model. There are three well known models for 
an isotropic linear viscoelastic liquid. In the Maxwell 
model a viscoelastic element is a combination of a purely 
elastic spring with a purely viscous dashpot. In the Bing- 
ham two-component model these elements are connected 
parallel to each other, and the Jeffreys element is a su- 
perposition of the two previous ones. A comprehensive 
review of the theoretical basics and the existing models 
was done by Bird in Refs. 00- An overview of the the- 
oretical approaches to the polymer dynamics is given by 
Doi and Edwards^. 

A ferrofluid with chains in a magnetic field is an 



anisotropic system. This makes it similar to a ne- 
matic liquid crystal (NLC), whose rheology can be de- 
scribed by a group of models usually referred as Ericksen- 
Leslie-Parody (ELP) models. The stress tensor in Refs. 
Ifl llfTITlTIT^ involves five independent viscosity coefficients, 
which are scalar functions of the density and the temper- 
ature. ELP models are widely used in the hydrodynamics 
of NLC. These models were employed in the analysis of 
a flow instability^ and in studying period ic pattern 
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In these papers as well as in Refs. Il5ll6l the influence 
of an external field was taken into account. Readers who 
are interested in a deeper insight into ELP models are re- 
ferred to the books of de Gennes^l and Chandrasekhai4&. 

Since the susceptibility of ferrofluids is much higher 
than that of liquid crystals, the effects caused by external 
field should be more pronounced. The theoretical analy- 
sis of the rheology of a colloidal suspension containing el- 
lipsoidal particles in a field was performed by Pokrovskij 
m Refs. Ilffctl The model gives a constitutive equation 
of the ELP kind, where the viscosity coefficients are ex- 
pressed in terms of the parameters of the suspension and 
the applied field. Electro- and magnetorheological flu- 
ids present a group of suspensions, whose behaviour is 
close to that of ferrofluids. The viscoelastic properties of 
electro- end magnetorheological fluids were studied ex- 
perimentally in Refs. l2lU22ll23l Particularly it was shown 
that external fields can influence the rheological response 
of the fluid by changing the complex shear modulus and 
yield stress. Zubarev and Iskakova™ used the results of 
Refs. ll9H20l to obtain a constitutive equation for a fer- 
rofluid with chains for the case of a weak flow of any 
kind. In the present paper the model suggested in Ref. 
13 is applied to the linear stability analysis of the Faraday 
instability. 

The Faraday instability denotes the parametric gen- 
eration of standing waves on the free surface of a fluid 
subjected to vertical vibrations. The study of this phe- 
nomenon dates back to the observations by Faraday^ in 
1831. The initially flat free surface of the fluid becomes 
unstable at a certain intensity of the vertical vibrations 
of the whole system. As a result of the instability, a 
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pattern of standing waves is formed at the fluid surface. 
The typical response is subharmonic, i.e., the wave fre- 
quency is half the frequency of the excitation. A har- 
monic response can be observed on a shallow fluid at low 
frequencies^,. Faraday waves allows one to investigate 
symmetry breaking phenomena in a spatially extended 
nonlinear system. Therefore they experience a renewed 
interest in recent years. Detailed experimental stud- 
ies of the various patterns on a viscous fluid have been 
performed 27 ' 28 Si22i2i, where a one-frequency as well as 
a two-frequency forcing were applied. Among the ob- 
served patters are parallel rolls 27 , hexagons 27 , a twelve- 
fold quasi-patternSi, triangles^, superlattices formed by 
small and large hexagons 29 , squares^&^iLiffli^I, and rhom- 
boid pattern^. 

The comprehensive linear stability analysis of the Fara- 
day instability on an arbitrarily deep layer of a viscous 
non-magnetic fluid has been performed by Kumar and 
Tuckerman 32 , This analysis was tested experimentally 33 
and an excellent agreement between the predicted and 
experimental data was found. In Refs. I26H341 the low fre- 
quency region is studied particularly. Bicritical points, 
where transitions from one type of response to others oc- 
cur, are predicted and experimentally confirmed— In 
Ref. 1351139 an analogy between the Faraday instability 
and a periodically driven version of the Rayleigh- Taylor 
instability is exploited. Based on that analogy in Ref. 
l35l a scaling law is suggested, which satisfactorily de- 
scribes the behaviour of the system in a wide range of 
parameters. S. Kumar— discusses the mechanism of the 
wave number selection in the Faraday instability on high- 
viscous fluids. 

In our previous paper—, the Faraday instability on a 
chain-free ferrofluid was studied. A nonmonotonic de- 
pendence of the stability threshold on the magnetic field 
is found at high frequencies of the vibrations. It was 
revealed that the magnetic field can be used to select 
the first unstable pattern of Faraday waves. In particu- 
lar, a rhombic pattern as a superposition of two different 
oblique rolls can occur. 

The Faraday instability of a viscoelastic non-magnetic 
liquid was studied experimentally in Ref. | 38l w here a 
harmonic response was detected. In Refs. 391.40l4ll the 
Maxwell model of viscoelastic liquid was used in the 
theoretical analysis. The authors observed pronounced 
changes in the neutral stability curves. Particularly, the 
tongues related to the harmonic response became abnor- 
mal. Such a tongue has no tip and all tongues of higher 
order are inside this abnormal tongue. 

The aim of the present paper is to investigate the role 
of the chains in the Faraday instability on a ferrofluid in 
a horizontal magnetic field. The formation of the chains 
leads to a dramatic increase of the magnetization relax- 
ation time, changes the effective viscosity of the suspen- 
sion, and increases the susceptibility of the ferrofluid. A 
number of model ferrofluids are investigated in a wide 
range of the parameters of the system to study the rela- 
tive importance of those effects. 



II. SYSTEM AND BASIC EQUATIONS 
A. Model 

A dielectric, viscous, and incompressible magnetic fluid 
with constant density p is considered, which contains par- 
ticles of the equal size. The strength of the dipole-dipolc 
interactions is characterized by the coupling constant 
e = /ioTO 2 /(167ri? 3 fcsT), which is the ratio of the en- 
ergy of interaction between two particles at the minimal 
separation with head-to-tail orientation of their magnetic 
moments to the thermal energy. Here m — 4Af s 7ri? 3 ) /3 
is the magnetic moment of a particle, R = Rq + S is 
the hydrodynamic radius of the particle, M s is the mag- 
netization of the magnetic material, Rq is the radius of 
the magnetic core of a particle, S is the thickness of the 
nonmagnetic layer, ks is the Boltzmann constant, and 
T is the temperature. The interaction of a particle with 
the applied field H is measured by the Langevin param- 
eter k — nomH J '(UbT). In the case of a magnetic fluid 
with low particle volume fraction ip, interactions between 
chains can be neglected. Assuming that the chains are 
straight and rigid, Zubarev and Iskakova determined the 
size distribution of chains, which minimizes the free en- 
ergy of such a magnetic fluid and is given bji24: 
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where v is the volume of a particle and y — mpexp{e). 
Each chain containing n particles is modelled by an uni- 
axial ellipsoid with semi-axes equal to nR and R. This 
keeps the solid phase volume density unchanged and al- 
lows one to use the results of PokrovskijiSi— The viscous 
stress tensor for a system of uniaxial ellipsoids consists 
of a symmetric and antisymmetric part 
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Both parts contain the geometric factors a n , f3 n , Xn, A„, 
p n , and £„, which are entirely determined by the as- 
pect ratio of an ellipsoid corresponding to a n-particle 
c j la j n i9 4 20 J 24 j 42 ^ j g v i scos ity of the carrier liquid, 

lik = (dkUi + diU k )/2 and Q ik = {d k Ui - diU k )/2 are 
the symmetric and antisymmetric parts of the tensor of 
velocity gradients, e and h = H/H are the unit vectors 
along the chain axis and the magnetic field, and the fol- 
lowing notations are used: 




n 

(• • • }n = / • • • &l) n (e)de, 



(•••)°= /•••e<(e)de. 



FIG. 1: A horizontally unbounded ferrofluid layer (2) is placed 
in a nonmagnetic container (3) with air (1) above. The system 
is subjected to a horizontal magnetic field Ho and harmonic 
vertical vibrations a(t). 



Here, ^(e) and ip n (&) are the equilibrium and the 
nonequilibrium angular distribution functions. The for- 
mer is known due to the classical Langevin model and 
the latter is the so lution of the Fokker-Planck equation 
(see Refs. Il9l20l24j) . The exact solution for the case of 
an arbitrary field is unknown. Therefore, Zubarev and 
Iskakova suggested an approximation for ^ n (e): 



^„(e) =V£(e) 



1 + a, 



(3) 



where a and b are a vector and a symmetric tensor. 
They can be found from the equations for the first (e k } n 
and the second (eie k ) n moments of ip n (e), which are de- 
rived from the Fokker-Planck equationi2*22iS4 with the 
accuracy up to the linear terms with respect to the ve- 
locity gradients. The moment equations for (e k ) n and 
(e-i£k)n involve two relaxation times, t% — l/(2D) and 
T2 = l/(6-D), where D ~ 1/r) is the coefficient of rota- 
tional diffusion. 

One can easily see, that the magnetization of the fer- 
rofluid is 
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R 3 
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If a and b are known one can relate the viscous stress 
tensor <j and the magnetization to the tensor of velocity 
gradients for the given value of magnetic field strength 
H. The total stress tensor T reads 



T lk = -\P + VoJ o MdH' + ^-H 2 j S tk + H t B k + o lkl 

(5) 

where B = /i (H + M), p is the pressure, and B the 
induction of the magnetic field. 



B. The system 

The above model is used to analyze the stability of the 
free surface of a ferrofluid in the following setup. The 
laterally infinite ferrofluid layer of arbitrary depth d is 
subjected to a homogeneous dc horizontal magnetic field 
Ho = (Ho, 0, 0) and harmonic vertical vibrations (Fig.^l. 
The plane 2 = coincides with the nondeformed surface 
of the ferrofluid. The fluid layer is bounded from below 
by the bottom of the nonmagnetic container and has a 
free surface described by £(f, x, y) with air above. 

Due to zero electrical conductivity of the fluid, the 
static form of the Maxwell equations is used for the mag- 
netic field in all three media. The fluid motion is gov- 
erned by the continuity equation and the conservation 
law of the linear momentum 



div u = 



— + (ugrad) m 



pg(t)6. 



(6a) 

i = 1,2,3. 
(6b) 



The vertical vibrations add a periodic term to the grav- 
ity acceleration go, i.e., a modulated value g(t) = 
(0,0,— go ~ a cos(wi)) appears in the equations of the 
fluid motion. Here a is the acceleration amplitude and to 
is the angular frequency of the vibrations. The governing 
set of equations has to be supplemented b y th e boundary 
conditions, which are the same as in Ref. 1371 



III. LINEAR STABILITY ANALYSIS 

Following the standard procedur e) 32 ! 43 , the governing 
equations and the boundary conditions have been lin- 
earized in the vicinity of the nonperturbed state 
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where Mq and pg are the magnetization and the pres- 
sure in the unperturbed state. In order to construct 
the linearized governing equations for the small perturba- 
tions, which encompass the magnetic field strength Hi, 
the pressure pi, and a nonzero velocity u of the fluid, it 
is necessary to expand all quantities in Taylor series. 



The above defined distribution functions g n and ip n 
are affected by the small perturbations. Since the forma- 
tion and the dissociation of chains are connected with the 
diffusion of particles in the suspension, they are rather 
slow processes. It can be estimated by the Schmidt num- 
ber Sc = r//(pDo) which relates the characteristic time 
for mass transport by flow to the characteristic time for 
mass transport by diffusion. For a typical fcrrofluid the 
Schmidt number is about 5 • 10 6 with 77 = 0.1 kg/(ms), 
p = 1020 kg/m 3 , and the Brownian diffusion coefficient 
Do = 2-10~ n m 2 /s — . Therefore we neglect the changes 
in g n caused by the perturbations. This implies that the 
size distribution of chains does not de pen d either on the 
spatial coordinates or on time. In Refs. 19 20. 2j| only spa- 
tially homogeneous systems are considered. An extension 
to the case of a spatially inhomogeneous system results 
in additional convective terms in the dynamic equations 
for ip n [Eq. (12) in Ref. |24|. Since these terms are of 
the form (ugrad) tpn,i, i-e., they are of second order with 
respect to the perturbations, they can be neglected in a 
linear stability analysis. As long as the typical length 
scale of a chain (~ 10 -8 . . . 10~ 7 m) is smaller than that 
of the spatial variation of the perturbations (typically 
~ 10~ 3 . . . 10~ 2 m), Zubarev's model can be applied in 
the present setup. Our results (see Sec. IIVII suggest that 
this condition is always fulfilled. 
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where A\, A2, ■ ■ ■ , A\ 3 are known functions of the ap- 
plied magnetic field, the frequency of the vibrations, as 
well as of the parameters of the ferrofluid, and the equa- 
tion of continuity in the form j xx + j vy + "f zz = has 
been used. The functions Aj are obtained by means of 
averaging (with function ip n (e)) over the orientations of 
chains and averaging over the sizes of the chains (with 
function g n ). 

Thus, the set of governing equation for the first order 
of the perturbations can be written: 

divB^=0, rotH^=0, i = 1,2,3, (7a) 
divu = 0, (7b) 

Here Bi is the perturbation of the induction of the mag- 
netic field. 

The linearized boundary conditions differ from those 
in Ref. 37 by the condition for continuation of the stress 
tensor across the free surface of the fluid 



The symmetry of the system possesses a number of re- 
strictions on a and b . The change of y to —y and z to 
— z should not change the components a x , b xx , b yy , and 
b zz , whereas this transformation leads to the change of a 
sign of the components a y , a z , b xy , and b xz . This com- 
pletes the set of the moment equations derived from the 
Fokker-Planck equationiSiSSiSi, and gives us the values of 
a and b . 



n j ( T S - T u}) ~ Si,s<rA±i = i = 1, 2, 3 at z = , 

(8) 

where = d xx +d yy . In contrast to the previous case 3 ^, 
the first two equations in (JSJ) are independent of each 
other. 

The stability of the flat surface with respect to stand- 
ing waves is analyzed by using the Floquet ansatz for the 
surface deformations and the z-component of the velocity 



£(t,x,y) = sin(kr)e 



(s-\-iauj)t 



(9a) 



After a and b are found, we have for the perturbations 
Mi of the magnetization and those of the total stress 
tensor T {2] : 



u z (t,x,y,z) = sin(kr)e( s+i ^ t ^ w n (z)e m ^ , (9b) 



n— — 00 



where k = (k x ,k y ,0) is the wave vector, s is the growth 
rate, and a is the parameter determining the type of the 
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response. For a = the response is harmonic whereas for 
a = 1/2 it is subharmonic. Expansions similar to (|9bl) are 
made for all other small perturbations and are inserted 
into the linearized governing equations 0. The functions 
of the vertical coordinate in the Floquet expansion are 
of the form w n {z) = C w e ±qz , where C w is a complex 
amplitude of the corresponding quantity. The condition 
of reality for £(t, x, y) leads to the equations^ 

£_„ = C a = , (10a) 
£-n = £_i, a = 1/2. (10b) 

Inserting the ansatz into the governing equations (7J|, 
we obtain a set of algebraic equations for the amplitudes 
of all perturbed quantities. The solvability condition for 
this set gives us four values of the modified wave vector 
q. Hence, the general solution of the set of governing 
equations in the ferrofluid contains eight arbitrary con- 
stants. Two more arbitrary constants are the amplitudes 
of the perturbations of the field in the air above and be- 
low the fluid after applying the boundary conditions at 
z — > ±oo. Ten of the boundary conditions [Eqs. (JSJ) and 
(3.2a)-(3.2f) in Ref. l37| allow one to express all the per- 
turbed quantities in terms of the coefficients £„ which 
satisfy the equation 

oo 

(W n £ n - o£ n _i - o£ n+ i) e [s+'(«+»Ht = , (11) 

n— — oo 

where W n are rather complicated functions of the applied 
field, the frequency of vibrations, the depth of the fluid 
layer, and of the parameters of the ferrofluid. Here the 
functions W n depend additionally on the coupling con- 
stant e in contrast to the corresponding functions in the 
previous paper—. 

Equation l|llfl has to be satisfied for all times which 
implies that each term of the sum equals to zero. Us- 
ing the relations between £ n with positive and negative 
numbers i|10|) , one gets the set of equations 

Wo£ - ad - a£i = 0, a = , (12a) 
^o£o - - <i = 0, a = 1/2 , (12b) 
W n £ n - <z£ n -i - a£n+i = 0, 71 = 1, •••oo. (12c) 

A cutoff at n = N (in the present work N = 100) leads to 
a self-consistent equation for the acceleration amplitude 

„45.46 



a = \F(a, k,u),H ,rj,a, p, K,e)\ , (13) 

where F is a complex function expressed in terms of con- 
tinued fractions. Equation l|13|) can be solved numerically 
and gives the dependence of a on k at fixed parameters. 
The critical values of the acceleration amplitude a c and 
the wave number k c correspond to the absolute minimum 
of the curve a(k) at zero growth rate (s =0). 



IV. RESULTS AND DISCUSSION 

Figure|21presents marginal stability curves for a viscous 
ferrofluid at high frequency for two cases. The dashed 
lines are calculated with neglect of the magnetization re- 
laxation time, a = and 6=0, i.e. the average orienta- 
tion of chains immediately follows the field perturbations. 
The solid lines depict the neutral stability curves for the 
system out of equilibrium. The dependence of the accel- 
eration amplitude on the wave number for s — divides 
the phase space into regions, where the surface of the 
ferrofluid is stable or unstable with respect to paramet- 
rically driven standing waves. The principal data, which 
can be extracted, are the critical acceleration amplitude 
(scaled with go), the critical wave number (scaled with 
the capillary wave number k a — y/pg/cr), and the num- 
ber of the tongue to which they belong. The number of 
a tongue I (from left to right) is the order of response: 
the basic wave frequency related to the l-th tongue is 
u>i = Ilo/2. The odd and even tongues are the regions, 
where either a subharmonic or a harmonic instability de- 
velops. That relation for the different instability types 
holds for Newtonian ferrofluids 3 £ but experiences signif- 
icant changes if the ferrofluid contains chains. 

Due to the finite relaxation time the tongues are de- 
formed. All tongues but one are now shaped by a lower 
boundary, a pronounced tip, and a upper boundary. For 
the chosen set of parameters in Fig. [21 the fourth tongue 
becomes exceptionally deformed, since it has no tip and 
no upper boundary in the way all other tongues have. 
It is caused by the fact that the self-consistent equation 
for the acceleration amplitude (|13f) has always a solu- 
tion as k goes to infinity. To note the difference to all 




FIG. 2: Neutral stability curves for the excitation frequency 
/ = 100 Hz and layer depth d = 5 mm. Dashed (solid) lines 
are calculated with (without) taking into account deviations 
of the magnetization from equilibrium. The parameters of 
the fluid are r\ — 0.1 kg/(ms), e = 5, n = 0.5, a — 0.0265 
N/m, p = 1020 kg/m 3 , R = 10.5 nm, Av = 614.5 m" 1 , and 
the particle volume fraction tp — 0.082. 
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FIG. 3: The second and the fourth tongues of the neutral 
stability curves for the viscosities of the carrier liquid rj = 
0.136 kg/(ms) (a), rj = 0.138 kg/(ms) (b), and rj = 0.14 
kg/(ms) (c). The remaining parameters are the same as in 
Fig.H 



other deformed tongues, we call it a abnormal tongue 
because this tongue lacks two essential features in com- 
parison with all other tongues^. The typical arrange- 
ment of tongues in Fig. [21 can be generalized as follows: 
it is always an even tongue, Z a b no = 2iV, N — 1, 2, • ■ • , 
which becomes abnormal (here / a bno = 4). All tongues 
with I > Z a bno form separated pairs of two overlapping 
tongues (here I = 5, 6 and I = 7, 8) inside the abnormal 
tongue. All tongues corresponding to a higher order har- 
monic response, I — l a bno+2L, L = 1, 2, • • • , are "islands" 
of stability with respect to a response of the system with 
the frequency w a bno- In the domains, where the subhar- 
monic tongues overlap with the regions of the harmonic 
type of instability (gray regions in Fig.0), both types of 
instability can occur. 

Figure [3] illustrates the scenario of the transition be- 
tween the neutral stability curves with Z a b n o = 4 and 
^abno = 2 with an increase of the viscosity of the car- 
rier liquid. At the point, where two tongues touch each 
other, the amplitudes for the corresponding terms in the 
Floquct ansatz © become equal. The number of the 
abnormal tongue / a bno depends on the mean relaxation 
time of a chain to = 1/ ((2D)}. The latter can be varied 
by changing the viscosity of the carrier liquid rj. In the 




T) [m/(kg s)] 

FIG. 4: Dependence of the product uj^bnoTa (see text) on the 
viscosity of the carrier liquid. The remaining parameters are 
e = 5, k = 0.3, a = 0.0265 N/m, p = 1020 kg/m 3 , R = 10.5 
nm and ip = 0.082. 



figure^] the product w a bnoTb is presented as a function of 
■q. As the relaxation time increases with the viscosity, the 
corresponding number ? a b n o decreases: from i a b n o = 60 
for r\ = 0.001 kg/(ms) to ; abn o = 4 for r\ = 0.1 kg/(ms). 
Therefore sequential changes of the number of the abnor- 
mal tongue occur. The corresponding value of w a bno T o 
shows a jump-like decrease at those points, where i a bno 

drops tO Z a bno - 2. 

Both the overlapping of tongues and the appearance 
of an abnormal one were also observed for the Maxwell 
viscoelastic liqui d 39 i 40 , and it seems that these features 
are signs of a viscoelastic behaviour. One has to note 
that if the angular distribution of chains is assumed to 
be in equilibrium (dashed lines in Fig. i.e. tq = 0, 
the neutral stability curves are standard tongues without 
any qualitative changes in comparison with the case of 
a Newtonian ferrofiuid (studied in Ref. f37t) . It shows 
that it is the combination of chains and their nonzero 
mean relaxation time which causes the new features in 
the neutral stability curve. 

The dependencies of the critical wave number and 
the critical acceleration amplitude on the excitation fre- 
quency f = lu/2tt are presented in Fig. |3J The curves 
denoted by A are calculated for a nonequilibrium mag- 
netization and the curves B present the case, where the 
chains follow the field perturbations immediately. The 
solid lines correspond to the ferrofiuid layer with depth 
al = 5 mm, and the dashed lines are calculated for an 
infinitely deep layer. It is seen that at high frequen- 
cies the critical acceleration in the system out of equi- 
librium is lower than that in a model system, where the 
relaxation time is neglected. It is clear that with a de- 
crease of the frequency the deviations of the magnetiza- 
tion from an equilibrium and their importance become 
smaller. Therefore, one can expect that curves A and B 
coincide at low frequencies. However, this is observed 
only in the case of the infinitely deep fluid layer. In 
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FIG. 5: Frequency dependencies of the critical wave number 
k c (a) and the critical acceleration amplitude a c (b) for d = 5 
mm (solid lines) and d = oo (dashed lines). The magnetiza- 
tion is assumed to be out of equilibrium for curves A and in 
equilibrium for curves B. The remaining parameters are the 
same as in Fig. [5] 

the case of d = 5 mm, the finiteness of the layer and 
consequently the influence of the viscous stresses in the 
bottom fluid layer become stronger in the case of the 
system out of equilibrium, and the critical acceleration 
increases more rapidly with a decrease of the frequency. 
With the further decrease of the frequency, transitions 
to higher order response occur^ 3 ^ 7 ., and the frequency 
of the Faraday waves remains high enough to leave the 
orientation of chains in nonequilibrium. 

Fig. presents the relative differences between the 
critical parameters fci free ^ and at? Tee ^ for a chain-free fer- 
rofluid and those for a ferrofluid with chains as functions 
of the Langevin parameter, i. e., of the dimensionless 
magnetic field. The magnetization, chord %c = M/H, 
and differential susceptibility \d — dM/ dH of the chain- 
free ferrofluid are chosen the same as in the ferrofluid 
with chains at each point. The magnetization is out of 
equilibrium in the ferrofluid with chains and is in equi- 
librium in the chain-free one. It is seen, that at / = 20 
Hz the formation of the chains decrease the critical wave 
number [curve A in Fig. Efa)]. Since the stresses at the 
bottom layer are essential at this frequency and d = 2 
mm, this leads to the increase of the critical acceleration 
amplitude [curve A in Fig.^fb)]. At the higher frequency, 
/ = 100 Hz, the changes of k c are non-monotonic. There 
is a range of the parameters, where the critical wave num- 
ber on the ferrofluid with chains is larger than that on 
the chain-free ferrofluid. The critical acceleration ampli- 
tude is decreased by the formation of chains [curve B in 
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FIG. 6: The relative differences between the critical wave 
numbers (a) and the critical acceleration amplitudes (b) as 
functions of the Langevin parameter. The quantities with the 
superscript "(free)" belong to the chain-free ferrofluid (g n = 
5 n ,i<p/v), the quantities without a superscript belong to the 
ferrofluid with chains [g n is given by Eq. Q]. Curve A: f = 
20 Hz, curve B: f = 100 Hz. The layer depth is d = 2 mm 
and the remaining parameters are the same as in Fig. [5] 



Fig. EP>)]. Note, that the susceptibility of the ferrofluid 
changes strongly with the increase of the field, therefore 
the relative importance of the influences of the magnetic 
field, the fluid microstructure, and the viscous stresses 
depends on the field. 

The above results were obtained for the particular case 
k||Ho. The typical size and orientation of chains depend 
on the applied field. Therefore, the threshold of the in- 
stability depends on the strength of the field and the 
angle 9 between k and Ho separately. It is not possible 
to introduce a single parameter like the effective field 37 , 
which would incorporate those two dependencies. This 
fact is illustrated in Fig. [7| The product H cos 9 was 
varied in two ways. For the curves indicated by A, the 
field strength was changed from zero to i?omax = 1-46 
kA/m with 9 = 0. In the second case (curves B), 9 was 
changed from zero to 7r/2 at the constant field i?omax- 
In contrast to the chain-free ferrofluid studied in Ref. 
l37l where both curves would be identical, there are now 
two different graphs. Whereas the critical wave num- 
bers have similar values for the two ways of variation, 
the critical acceleration shows greater differences: for in- 
stance for H cos9/H 0max — 0.5 the difference in k c is 
0.4 % against 2.5 % for a c . Thus it makes a notable 
difference for the stability of the free surface whether 
the system is subjected to (9 = 0,H = H 0max /2) or 
to (9 = 60°, H = H 0max ). This fact suggests a simple 
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FIG. 7: The dependencies of the critical wave number (a) 
and the critical acceleration amplitude on the projection of 
Ho on the direction of k for / = 100 Hz and d = 1 mm. 
The curve A is obtained by varying the field strength with 
9 — 0, whereas the curve B was calculated by varying 6 with 
Ho ~ -ffomax = 1-46 kA/m. Remaining parameters are the 
same as in Fig. [5] 



method to test the presence of chains in a sample. If 
the fcrrofluid is chain-free both ways of excitation lead 
to the same value of a c , if it contains chains one gets two 
different values. 

The dependencies are less pronounced than in the pre- 
vious case, where -ffomax = 26.6 kA/m2I, due to the fact 
that the range of rather weak fields is studied here. Such 
a choice of fields ensures that (n) < e and that the head- 
to-tail orientation of the moments in a chain is still prefer- 
able which is needed to apply the model. 

The nonmonotonic dependence of a c (Ho) for 8 = 
[curve A in Fig.[7fb)] caused by the joint action of the two 
mechanisms of the viscous damping: i) the dissipation in 
the bulk fluid and ii) the viscous stresses in the bottom 
fluid layer. The first mechanism is dominant for large 
wave numbers and causes a decrease of a c . The second 
becomes important with the decrease of the wave num- 
ber and results in an increase of a c with field strength. 
A detailed discussion of this nonmonotonic behaviour of 
a c (H) is given in Ref. l3?l 

There are three effects caused by the formation of 
chains. Since a chain has a magnetic moment larger than 
that of a single particle, the formation of chains leads to 
the increase of the magnetization of the system at the 
given field in comparison with a homogeneous ferrofluid. 
This strengthens the field influence on the fluid motion. 
The effect can be formally taken into account by adjust- 



ing the magnetization and susceptibility of the chain-free 
fcrrofluid in such a way as it is done in Fig. The two 
other effects are more interesting. The chains increase 
the magnetization relaxation time and change the effec- 
tive viscosity of the suspension both in a steady and in a 
periodical flov,r 19 i 20 i 24 . 

In order to investigate the role of these two effects, 
a number of model ferrofluids are considered. The first 
sample is a ferrofluid with chains. It is characterized by 
the mean relaxation time of a chain tq and unperturbed 
values of the magnetization Mo as well as the chord Xc 
and the differential susceptibility \d- The second model 
ferrofluid contains as well chains, but in this case the re- 
laxation time is neglected [VVi(e) = "0n( e )]- The third 
and the fourth ferrofluid is chain-free (g n = S n ,iip/v) 
with the same unperturbed values of Mo, Xc, and \d- 
The magnetization relaxation time for the third ferrofluid 
is equal to To, and for the forth sample the relaxation 
time is neglected, i.e. the last sample is the ferrofluid, 
which was studied in Ref. IU The predictions of the 
Zubarev model for the case of a chain-free fcrrofluid with 
a finite relaxation time were compared with the results 
of a model using the effective field theory developed by 
Martsenyuk, Raikher, and Shliomis in In the 

both cases, the Langevin parameter K a dj involved in the 
model was adjusted to obey mip/vL(n a dj) = Mo, where 
L(k) = coth(K) — 1/k denotes the Langevin function. 
A fairly good agreement between the threshold param- 
eters predicted by both models is observed. In Fig- 
ure [S] the coupling constant e is varied by changing the 
temperature. The intensity of the applied magnetic field 
Hq is changed in such a way that the Langevin parameter 
is equal to k = 0.5. 

The comparison of the presented dependencies reveals 
that the finiteness of the magnetization relaxation time 
leads to an increase of the critical wave number [compare 
the pairs of the curves A - C and B - D in Fig. Efa)]- 
At the same time, the changes in the slow-flow rheologi- 
cal properties of the suspension caused by the formation 
of the chains decrease strongly the critical wave number 
(compare curves C - D and A - B). The net effect of 
those two influences gives the behaviour that is presented 
by the curves B in Fig. |S] 

The critical acceleration amplitude for a ferrofluid 
with nonequilibrium chains [curve B in Fig. IHJb)] has 
the lowest value among all the tested model ferroflu- 
ids. This can be interpreted as a sign of an elastic be- 
haviour, which is typical for viscoelastic liquids at high 
frequencies. This fact is confirmed in Fig. where the 
real and the imaginary parts of the effective viscosity 
VcS = [fife / ip'likfi) ~ °f a ferrofluid are presented as 
functions of the frequency of surface waves /; = u>i/(2tt) 
with 1 = 1, The case Hq = is considered, where the 
system is isotropic and the influence of the chains on the 
rheological properties of the suspension can be character- 
ized by the single complex parameter rj e f^. The real part 
of the effective viscosity determines the dissipation of the 
energy in the system, and a nonzero imaginary part leads 
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FIG. 8: The dependencies of the critical wave number (a) 
and the critical acceleration amplitude for the infinitely deep 
layer of the ferrofiuid consisting of A: single particles with 
the Brownian relaxation time tb = to; B: chains with a 
orientations out of equilibrium [Eq. J3J]; C: single particles 
with zero relaxation time; D: chains with the equilibrium 
orientations [ip n (e) — ip^(e)]. The remaining parameters are 
the same as before. 



0.001 1 — — 1 

1 10 100 

/,,Hz 

FIG. 9: The dependencies of the real (a) and the imaginary 
(b) part of the effective viscosity on the wave frequency. The 
dashed line denotes the effective viscosity, which is calculated 
with neglect of the relaxation time for e = 5. The solid lines 
denote the effective viscosity calculated using the size distri- 
bution of chains l(T)l. and the dotted curves present the real 
and imaginary parts of the effective viscosity for model fer- 
rofluids containing chains of n particles only. Ho = 0, the 
remaining parameters are the same as before. 



to a phase shift between the tensor of velocity gradients 
and the stress tensor, i.e. to a viscoelastic behaviour. It 
is seen that the imaginary part of the effective viscosity 
is nonzero and has a maximum at a frequency of /; ~ 5 
Hz [solid line in Fig. EJb)], where therefore a maximal 
phase shift (viscoelastic behaviour) can be expected. At 
the same time, the real part of the effective viscosity de- 
creases [solid line in Fig.|^|[a)]. The effective viscosity for 
a stationary flow (with neglect of the relaxation time) is 
plotted as dashed line. It is seen that the deviations of 
the orientational distribution from equilibrium result in 
a decrease of the real part of the effective viscosity. This 
implies that the energy is dissipated weaker, and that 
consequently a c is lower than in the case of zero relax- 
ation time. The comparison of the solid curve with the 
dotted ones in Fig. |5Jb) shows that the maximum of the 
imaginary part of the effective viscosity is caused by the 
long chains. 

A comparison of the contributions of chains of different 
size is done in Fig. ^| for three values of the wave fre- 
quency. The size distribution of chains Q is presented 
with filled circles in part (a) of that figure. It is seen 
that at low frequencies (/ = 1 Hz) the viscoelastic effect 
is mainly caused by very long chains with small values of 
g n [see maximum of —Im(r/ C g n ) at n = 16 in Fig. llUf b)]. 
For / = 5 Hz the maximum of — Im(r? e fltn) is shifted to 
long chains (n = 11) and the value of the maximum in- 
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FIG. 10: Contributions of chains containing n particles in the 
real (a) and in the imaginary part (b) of the effective viscosity 
at fi = 1 Hz (o), 5 Hz (□), and 100 Hz (0). The filled circles 
present the size distribution of chains (0. The remaining 
parameters are the same as in Fig. [5] 
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creases. Together with the larger values of g n for those 
long chains, the total viscoelastic effect is larger than for 
/ = 1 Hz [compare the area under the o- and D-curve in 
Fig. llOf b)]. For high frequencies as / = 100 Hz the max- 
imum of — Im^cffn) is at short chains (n = 5) and the 
value of the maximum decreases. This decrease cannot 
be compensated by the even larger values of g n for those 
short chains. Thus the total viscoelastic effect decreases 
again: the area under the ^-curve is smaller than the 
area under the D-curve. 

In summary, by increasing the frequency from 1 to 100 
Hz the viscoelastic effect passes at maximum at / = 5 
Hz [solid curve in Fig.Efb)]. The maximal contribution 
to — Im(ry e fj n ) comes for / < 100 Hz from chains off the 
maximum of g n . Comparing the area under the o-, □-, 
and ^-curve in Fig. HOf ah one can see that the area is 
continuously shrinking with /. It results in the solid line 
in Fig. Eta) describing the weaker dissipation of energy 
in the system. 

The origin of the viscoelasticity (i. e., the dependence 
of the relation between the tensors of velocity gradients 
and the stress tensor on the history) in a ferrofluid with 
chains can be explain in the following way. When the 
orientation of a chain is changed from the equilibrium 
by the velocity gradients, a restraining moment is gener- 
ated by the applied field. Thus the network of chains is 
an "elastic" element of the system, which interacts with 
the viscous carrier liquid. If the relaxation time of a 
chain is much shorter than the period of the vibrations, 
the orientation of chains and consequently the rheologi- 
cal properties of the suspension do not depend on history, 
and there is no viscoelastic properties. This fact becomes 
more obvious, if one notes that the relaxation time of a 
chain is proportional to the carrier liquid viscosity. Thus 
the negligible relaxation time of a chain implies a weak 
coupling between the chain and the carrier. 

To estimate the importance of the relaxation processes 
in a ferrofluid with chains, the Deborah number De is 
calculated for the subharmonic response. The Deborah 
number is the product of the angular wave frequency 
with the mean relaxation time of a chain, De = loiTq. 
In the case of polymer solutions, the most pronounced 
viscoelastic behaviour should be expected in the range of 
De w 1, whereas at De > 1 a recovery of the Newtonian 
behaviour was observed^. In the case of ferrofluids with 
chains, the viscoelastic behaviour is present over a much 
wider De range. The frequency range of Fig. E] gives 
0.033 < De < 3.3 for which a viscoelastic behaviour 
is present since — Im(r/ ff) ^ 0. That fact is supported 
by the range of appearance for the abnormal tongues 
considered to be a characteristic feature of viscoelastic 
behaviour. From Fig. 0] one can read that viscoelastic 
properties can be expected even for De beyond 10. 

Obviously the viscoelasticity of a ferrofluid with chains 
differs from that of a polymer solution. The reason for 
this difference is the following. The relaxation time of 
a long polymer chain involved in the polymer dynam- 
ics is the conformational relaxation time. As soon as 



the period of excitation becomes much shorter than the 
conformational relaxation time, the polymer chain does 
not participate in the movement of the solution. There- 
fore a motion of the Newtonian solvent through the rigid 
obstacles of polymer chains is observed. In the case of 
the rather short chains formed by magnetic particles, the 
relaxation time is the orientational relaxation time of a 
chaini^ in the viscous carrier liquid. Since it is propor- 
tional to the viscosity of the carrier liquid, the high val- 
ues of t and De, respectively, mean a strong coupling 
between the motion of a chain and the carrier liquid, 
which ensures the viscoelastic behaviour. Thus ferroflu- 
ids widen significantly the range in which viscoelastic fea- 
tures can be observed in contrast to polymer solutions. 

In a real ferrofluid, any variation of the system param- 
eters not only varies the Deborah number but also af- 
fects the critical parameters either directly or indirectly 
by changing other properties as permeability or viscosity. 
In order to investigate the role of the deviations of the 
chains from equilibrium in greater detail, a model fer- 
rofluid with a tunable coefficient of rotational diffusion 
for all chains D tunc ^ n — D n D scai \ e is studied. Via that 
coefficient the mean relaxation time To, i.e. the Deborah 
number, can be varied. 

-Dscaic is a phenomenological parameter, which is re- 
lated neither to the chain size nor to the fluid properties. 
Using this phenomenological parameter, a smooth tran- 
sition can be realized from a model ferrofluid with chains 



i-6 - ^<^yv 

(a) 

1.2 





26 








© 
bo 


24 




\c 


(b) - 




22 










20 









0.001 0.01 0.1 1 10 

De 



FIG. 11: Dependency of the critical wave number (a) and the 
critical acceleration amplitude (b) on the Deborah number 
for a model ferrofluid with the tunable relaxation time of a 
chain. The coupling constant and the Langevin parameter 
are: Curves A: e = 5.5, k — 0.39, B: e — 5, « = 0.5, and C: 
e — 4.5, k — 0.63. The remaining parameters are the same as 
before. 
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almost in equilibrium (-D S caio ^ 1 gives De <C 1) to the 
real ferrofluid (-D sca ie = 1); an d further to a model fer- 
rofluid with chains far from the equilibrium (D scll \ c <C 1 
gives De 3> 1). 

It is seen in Fig. ^] that even small deviations from 
equilibrium are significant as the variation of k c and a c 
in the range from De ~ 0.001 to De ~ 0.01 show. In the 
case of a system with long chains (e = 5.5), a pronounced 
minimum in the dependence of a c (De) and a maximum 
in that of k c (De) are observed. As e and the mean chain 
length decrease, these extrema are shifted towards higher 
values of De and become less pronounced and disappear 
for e = 4.5. Beyond De ~ O(l) one observes a saturation 
of k c and a c close to the corresponding extremal values. 
That means that beyond a certain deviation from equi- 
librium, viscoelastic effects can neither be enhanced nor 
diminished. The latter confirms that there is no recovery 
to a Newtonian behaviour at high De numbers. 

V. CONCLUSION 

The rheological properties of a ferrofluid caused by the 
formation of chains is investigated by the linear analy- 
sis of the Faraday instability in a horizontal magnetic 
field. A horizontally unbounded ferrofluid layer of a fi- 
nite depth has been considered. The dependencies of the 
critical acceleration amplitude a c and the critical wave 
vector on the excitation frequency /, the magnetic field 
Ho, and the dipolar coupling constant e have been ob- 
tained for different depths of the layer in a wide range 
of fluid viscosities. A viscoelastic behaviour is predicted, 
which is indicated by the qualitative changes in the neu- 
tral stability curve: an abnormal tongue appears which 
has no tip and upper boundary in contrast to all other 
tongues (Fig. [5J. Beside its existence it is revealed how 
this abnormal tongues is formed by the merging of other 
tongues as the viscosity of the carrier liquid is changed 



(Fig.©. 

The threshold of the instability depends on the applied 
field and on the angle 9 between the wave vector and 
the applied field separately (Fig. Q) . Therefore one can 
easily test a ferrofluid whether it contains chains or not by 
choosing different combinations of 9 and Hq but the same 
value for Hq cos 9. If the ferrofluid contains chains one 
get different thresholds whereas the surface of a chain- 
free ferrofluid is destabilized at a unique value of a c . 

A ferrofluid with chains whose orientation is out of 
equilibrium compared with other model ferrofluids (Figs. 
El El has the lowest threshold of all. This is caused by 
a decrease of the real part of the effective viscosity [Fig. 
EJa)] which corresponds to a weaker dissipation of energy 
and thus to a lower value of a c . The nonzero imaginary 
part of the effective viscosity [Fig.^b)] indicates also the 
viscoelastic behaviour of the studied ferrofluid. By ana- 
lyzing the dependence of the critical parameters on the 
Deborah number, it is found that a viscoelastic behaviour 
occurs over a much wider range than in viscoelastic poly- 
mer solutions. Particularly no recovery to a Newtonian 
behaviour at high Deborah numbers is found. 

The present model is restricted to the case of a 
monodisperse ferrofluid. An account of the polydisper- 
sity will be necessary to compare the predictions with 
future experiments using real ferrofluids. The dipolar 
interactions between chains and their flexibility will in- 
fluence the properties of a real ferrofluid. Nevertheless, 
the above-discussed effects should be detectable, at least 
qualitatively, in a real experiment. 
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